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Abstract 



In culture migrating and interacting amoeboid cells can form ne- 
matic liquid crystal phases. A polar nematic liquid crystal is formed 
if the interaction has a polar symmetry One type of white blood cells 
(granulocytes) form clusters where the cells are oriented towards the 
center. The core of such an orientational defect (disclination) is ei- 
ther a granulocyte forced to be in an isotropic state or another cell 
type like a monocyte. An apolar nematic liquid crystal is formed 
if the interaction has an apolar symmetry. Different cell types like 
human melanocytes (=pigment cells of the skin), human fibroblasts 
(=connective tissue cells), human osteoblasts (=bone cells), human 
adipocytes (= fat cells) etc., form an apolar nematic liquid crystal. 
The orientational elastic energy is derived and the orientational de- 
fects (disclination) of nematic liquid crystals are investigated. The 
existence of half-numbered disclinations show that the nematic phase 
has an apolar symmetry. The density- and order parameter depen- 
dence of the orientational elastic constants and their absolute values 
are estimated. From the defect structure, one finds that the splay elas- 
tic constant is smaller than the bend elastic constant (melanocytes). 
The core of a disclination is either a cell free space or occupied by non 
oriented cells (isotropic phase) or occupied by a cell with a different 
symmetry or occupied by another cell type. 
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Introduction 



Investigations of the structure and physical properties of liquid crystals have 
greatly increased in the last decades and are now an integral part of solid- 
state physics. Certain organic materials do not show a single transition from 
solid to liquid, but rather a cascade of transitions involving new phases. 
The mechanical properties and the symmetry properties of these phases are 
intermediate between those of a liquid and those of a crystal 

Depending upon the nature of the building blocks and upon external 
parameters (temperature, solvents, etc.) a wide variety of phenomena and 
transitions amongst liquid crystals are observed. To generate a liquid crys- 
tal one must use anisotropic objects, like elongated molecules. Several ways 
are known to achieve this: with small molecules; with long helical rods that 
either occur in nature or can be made artificially; with polymers; with more 
complex structure units like the capside of a virus @ that are associated 
structures of molecules and ions; likewise, with even more complex structure 
units, like amoeboid cells which are complex structures far from thermo- 
dynamic equilibrium and act as a machine |3], f|]. Here, the orientational 
elastic properties of liquid crystal formed by migrating and living cells are 
investigated. 

Thermotropic nematic liquid crystals can be formed by small elongated 
molecules like p-azoxyanisole (PA A) or N-(p-methoxybenzylidene)-p-butylaniline 
(MBBA) which can be regarded as rigid rods (~ 2 nm length and ~ 0.5 nm 
width). At low temperature one observes a crystalline phase. Heating up the 
sample, a transition to a nematic phase is observed. With further heating 
the transition to an isotropic liquid is observed. 

Lyotropic nematic liquid crystal can be formed by long elongated molecules 
in suitable solvents. The anisotropic building blocks can have a rod-like con- 
formation like (i) synthetic polypeptides (~ 30nm length and ~ 2 nm width) 
and (ii) tobacco mosaic virus (~ 300 — 1000 nm length and ~ 20 nm width). 
At low concentration one observes an isotropic liquid but a transition to a 
nematic phase is observed by increasing the concentration of the elongated 
building blocks. 

The fundamental property of a nematic liquid crystal, which makes it 
different from an isotropic liquid and similar to a solid, is the presence of 
an orientational degree of freedom which is characterized by a macroscopic 
spatial ordering of the long axis of molecules. The statistical mechanics for 
predicting the macroscopic spatial ordering is difficult but can be approxi- 



3 



mated by applying mean-field calculations, (i) The Onsager approach starts 
with anisotropic steric repulsion. The free energy is minimum for oriented 
rods (width << length): they show an orientational order of a nematic phase 
at high rod density. But, at low rod density, the suspension is an isotropic liq- 
uid, as expected, (ii) The Maier-Saupe theory approaches the intermolecular 
interactions by a mean-field which is proportional to the order parameter of 
the nematic phase. The minimum of the free enthalpy yields the equilibrium 
state. At high temperature an isotropic liquid without orientational order is 
predicted, but a nematic liquid crystal with orientational order can occur at 
low temperatures. 

The above mentioned nematic liquid crystals are described by equilibrium 
thermodynamic. This holds true even for tobacco mosaic virus since only the 
anisotropic shape of the viruses is considered. The intracellular metabolism 
can be neglected since the virus does not have the ability of self-motion. 
Thus, the thermal motion is the stochastic source for the molecular dynam- 
ics. However, liquid crystal phases can also be formed if the stochastic source 
is created by energy consuming processes. A typical example are migrating 
biological cells like white blood cells. In the absence of a directing extra- 
cellular signal, the single cells perform a random walk due to the existence 
of stochastic processes in the intracellular signal transduction chain. The 
migrating cells can form liquid crystal phases if each cell transmits extracel- 
lular signals that are received by other cells. The symmetry of the cell-cell 
interaction is important for the formed type of liquid crystal phase: 

(i) A polar nematic liquid crystal is expected if the cell-cell interaction is 
polar. This means that one cell transmits a signal which attracts another cell. 
The cellular response is directed migration. For example: (a) Single slime 
mould cells migrate as amoeboid cells on a surface and search for bacteria as 
food. In the case of low bacteria concentration, the intracellular metabolism 
of the slime mould cells alters in such a way that cyclic AMP is emitted as 
an extracellular signal. Other slime mould cells can receive this chemical 
signal and react accordingly to the received signal with directed migration. 
At high enough cell density, the cells form a cluster where each cell tries to 
move towards the center of the cluster |J. (b) A similar process is observed 
with migrating white blood cells ||. If human leukocytes (granulocytes) 
are expose to blood plasma they then migrate as single amoeboid cells on a 
surface as shown in Fig. |l|a. But if the calcium concentration is lowered, e.g. 
by adding the chelator EDTA (ethylene diamine tetraacetic acid, pKc a = 
10.7) or EGTA (ethyleneglycol-(aminoethyl ether) tetraacedic acid, pK Ca = 
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11) to the plasma, then the migrating granulocytes attract each other and 
form a cluster, as shown in Fig. |]b. The emitted molecules (chemical signal) 
is still unknown. The process is reversible since the clusters disappear if 
the calcium concentration is increased. In both examples, the cells in a 
cluster form a polar nematic liquid crystal because (i) there is order in the 
cell orientation, (cosip), induced by the directed migration and (ii) no order 
in the center of mass of the cells. The angle, <p, is the difference between 
the actual migration direction and the direction towards the center of the 
cluster. The cells in the cluster of Fig. |I|b have a polar order of 0.82. A 
singularity in the cell orientation is expected in the center of the cluster. A 
common observation is that the cell in the center has a spherical shape and 
the surrounding cells have a polar shape. Human monocytes (another type of 
leukocytes) are immobile on a glass surface. If a mixture of granulocytes and 
monocytes is investigated, the immobile monocytes act as nucleation center 
as shown in Fig. The polar nematic liquid crystal formed by living and 
migrating cells is a system far from thermodynamic equilibrium. 

One can ask which temperature corresponds to the random movement 
of the cells? This question can be answered by using the Einstein relation 
which connects the random movement of an inert particle described by the 
diffusion coefficient D, the mobility R and the thermal energy ksT: DR = 
fcflT '. The mobility of a spherical particle with radius r in a viscous medium 
with viscosity 77 is given by Stokes: R = 6tn]r. The random movement of 
a migration granulocyte can be quantified by measuring the mean-squared 
displacement as a function of time. The diffusion coefficient, D, is one fitting 
parameter (D = 200 — 1000 fim 2 /min). The calculated temperature for an 
inert particle which has the value of the diffusion coefficient of the migrating 
cells, is very high w 10 5 K (r = 10/xm, rj = 0.01P) ||. 

(ii) An apolar nematic liquid crystal is expected if the cell-cell interaction 
is apolar. This means that one cell can influence the orientation of another 
cell. The cellular response is such that the cells like to be parallel. For ex- 
ample, human melanocytes are elongated cells which distribute with their 
dendrites the pigment melanin in the skin. Although the involved chemical 
signal is unknown there exists an attractive cell-cell interaction because the 
cell body of one cell attracts another cell body. There also exists a repulsive 
cell-cell interaction since the dendrites avoid contact with each other. The 
elongated cells try to be parallel oriented to each other and form at high cell 
density a nematic liquid crystal. At low cell density the cells form cluster 
with oriented cells as shown in Fig. |2|a, however, with increasing cell con- 
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centration the cells form a nematic liquid crystal. Fig. shows a nematic 
liquid crystal close to the nematic-isotropic transition. This observation is 
not restricted to melanocytes. Other cell types like fibroblasts, osteoblasts, 
adipocytes etc. orient in the same way so that, in the mean, the cells are 
parallel to each other. The cells form an apolar nematic liquid crystal be- 
cause there is no order in the center of mass of the cells but order in the 
cell orientation. The apolar symmetry of this phase will be demonstrated by 
investigating orientational defects as shown below. A few words to the dy- 
namic behavior: Time lapse movies of melanocytes show that the cell body 
performs rhythmic movements and, in addition, each cell extend and retract 
periodically its opposing dendrites. The dendrites find their orientation with 
minimum interaction with other cells J7J. 

A uniform oriented nematic liquid crystal is formed if a small extracellular 
guiding signal is applied like parallel scratches in the surface (Fig. |2|). With 
this method oriented single nematic liquid crystals were produced which had 
a size of several cm 2 . An isolated single cell is not oriented by such a weak 
extracellular guiding signal The uniform oriented nematic liquid crystals 
which are used in the liquid crystal display technology, are made in a very 
similar way M. In the absence of an orienting extracellular guiding signal the 
nematic phase contains numerous orientational point defects (disclinations). 
An example is given in Fig. |3|. 

Summarizing this part, nematic liquid crystal phases can be formed ei- 
ther by anisotropic interacting building blocks at thermal equilibrium or by 
anisotropic interacting building blocks far from thermodynamic equilibrium. 

Orientational Elastic Energy 

An apolar nematic liquid crystal is an uniaxial system ( ±n) where the 
molecules are on average aligned along the director ±n. In the case of an 
ideal nematic liquid crystal, the direction of the director ft is the same all 
over the sample. However, this ideal conformation will not be compatible 
with the constrains that are imposed by limiting boundaries. There will be 
some deformation of the alignment. The distance £, over which significant 
variations of orientation occur, is of most interest, and is much larger than 
the dimensions, a, of the building blocks. Thus, the deformations may be 
described by a continuum theory disregarding the details of the structure 
on the scale of the building blocks (a /I « 1). To construct such a theory 



6 



one possible starting point would be the elastic distortion energy density F d 
(H U). The spatial derivatives of n{f) are used to construct F d . They form 
a tensor of rank two c^ng. At any point we introduce a local system of 
Cartesian co-ordinates, x, y, z with y parallel to n at the origin, x chosen 
perpendicular to y (x-y plane are the plane where the cells migrate) and 
z parallel to the normal of the plane, x, y, z form a right-handed system. 
Referred to these axis, the two components of curvature, at this point, are the 
splay and bend deformation s = dn x /dx and b = dn x /dy, respectively (see 
Fig. £|). The distortion energy of a liquid crystal specimen in a particular 
configuration, relative to its energy in the state of uniform orientation, is 
expressible at the area integral of distortion energy density g which is a 
quadratic function of two differential coefficients which measure the curvature 
in two dimensions. 

g = k lS + k 3 b + T^his 2 + -ki 3 sb + -k 33 b 2 H (1) 

The spontaneous splay coefficient, ki is zero if the axis perpendicular to ft is 
a symmetry axis (•<—=—>•)■ The spontaneous bend coefficient, k 3 , is zero if the 
axis parallel to n is a symmetry axis (f=!)- The splay-bend coefficient, k 13 , 
is only unequal to zero for — and TtH- &11 an d ^33 are the orientational 
elastic coefficient for splay and bend deformation. These coefficients are also 
known as Frank elastic constants. 

The distortion energy can be written in the coordinate free notation 

11 1 
g = k l Vn + k 3 (curln) z + -k ll (Vn) 2 + -ki 3 (^ (2) 

Orientational Elastic Constants 
Splay and Bend Constants 

The orientational elastic constant of a nematic liquid crystal (bulk) can be 
estimated by purely dimensional arguments: one expects that the elastic 
constants are to be in the order of U / a, where U is a typical interaction energy 
between the molecules while a is a molecular dimension. A typical value of 
the molecular interaction energy is ~ 8 kJ/mol or 1.3 x 10~ 20 J/molecule 
which corresponds to a temperature of ~ 2000 K. A typical liquid crystal 
molecule has a length of about 1.5 nm. It leads to a value of ~ 8 x 10~ 12 
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N which is in accordance with the experiments The orientational elastic 
constant of a membrane can be estimated by taking this value times the 
thickness, d, of the membrane. It yields 3 x 10~ 20 J which is again verified 
by experiments [[UJ . 



The orientational elastic constant of a nematic liquid crystal formed by 
living cells can be estimated in a similar way: In the case of migrating gran- 
ulocytes the quasi temperature is ~ 10 5 K. A condensed phase can only be 
formed if the interaction potential, U, is in this order or even larger. This 
leads to an interaction energy of ~ 5 x 10 -19 J/cell. The bulk elastic con- 
stant is obtained by dividing this interaction energy by the length of a cell 

20pm). It yields 2.5 x 10~ 14 N. The two dimensional elastic constant is 
obtained by multiplying this value with the cell thickness (~ 20pm). The 
estimated orientational elastic constant of a quasi two dimensional nematic 
liquid crystal is then « 5 x 10~ 19 J. 

The dependence of the orientational elastic coefficients on the apolar order 
parameter, S, and the cell density, p, can be estimated in the following 
way: The orientation angle of amoeboid cells is controlled by an automatic 



controller which compares the actual angle with the desired one ||12|| . The 
cellular response is in such a way as to minimize the deviation between the 
actual angle and the desired one. The desired angle of one cell is to be parallel 
to another cell. In the case of a nematic liquid crystal the information of the 
desired angle is transferred to one cell by an extracellular guiding field which 
can be chemical, electrical or sterical. The extracellular guiding field, E 2 , of 
a nematic liquid crystal can be approached by the nematic mean field where 
one cell is considered in the field produced by all the other cells. The strength 
of the mean field is assumed to be proportional to (i) the cell density, p, and 
(ii) the mean apolar order parameter, S(= (cos20)), of the surrounding cells. 

E 2 = a 2 pS (3) 

These assumptions are verified by experiments ||. The extracellular guiding 
field is calibrated in cellular units by a 2 . 

The cellular automatic controller acts in such a way as to minimize the 
deviation between the actual angle and the desired one. In the absence 
of extracellular guiding field, the cells alter their angle of orientation in a 
random fashion. The response of the cellular automatic controller is hence 
not simply described by the extracellular guiding field but also by stochastic 



processes in the cellular signal transduction chain [12]. The rate equation for 
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the orientation angle is a stochastic differential equation (Langevin equation): 

dQ 

— = -k 2 E 2 sin2G + r 2 (t) (4) 

The orientation angle, 0, is measured in respect to the director n of the ne- 
matic liquid crystal. The cellular signal transformer is described by E 2 sin2Q 
and the cellular reaction unit by the coefficient k 2 . The stochastic part of 
the machine equation, T(t), is approximated by a white noise of the strength 
q 2 ((r) = 0, (T(t) r(t')) = q 2 5(t — t')) . The Langevin equation can be trans- 
formed into a Fokker-Planck equation which describes the temporal varia- 
tions of the angle distribution function ]13|, ||. The predicted steady state 
angle distribution function is 

/ 2 (e) = / 20 e yc - 2e (5) 

with the generating function V 

2 k 2 E 2 2k 2 a 2 p S 
V = = = A 2 p S (6) 

92 92 

f 2 o is determined by the normalization (/ f(Q) dQ = 1). It is important to 
realize that the steady state angle distribution function of living cells (system 
far from thermodynamic equilibrium) have the same mathematical structure 
as the Boltzmann distribution which describes the fluctuations of a system at 
thermal equilibrium. The predicted angle distribution function, /2(@), of a 
nematic liquid crystal formed by migrating and orienting cells as well as the 
density- and order parameter dependence of the extracellular guiding field 
are verified by experiments Q. The unknown coefficient, A 2 (= 2 k 2 a 2 / q 2 ), 
can be determined by measuring the angle distribution function at different 
cell densities. For melanocytes one gets l/A 2 = 55 cells/mm 2 . Form the 
self-consistency for the order parameter one gets the transition from the 
nematic to the the isotropic phase ||]: A 2 p = 2 with cell density, po, at the 
transition. This threshold cell density, po, can be estimated by considering 
random oriented cells without steric contact. One gets a threshold cell density 
of ~ 100 cells/mm 2 for »s 100/im long cells. This value is in accordance with 
Po (= HO cells/mm 2 ) determined directly from A 2 . 

The interaction potential, U, of a single cell with its environment is the 
generating function, V, times cos2G 

U = A 2 pS cos26 = 2— S cos20 (7) 

Po 
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This interaction potential of a single cell in the environment of a nematic 
liquid crystal will be used to predict the order parameter dependence as 
well as the density dependence of the orientational elastic coefficients. The 
calculation is made in analogy to the elastic coefficients of nematic liquid 



crystal formed by molecules 14 



The angle alters by the distortion angle a if one moves in space (y- 
direction). One obtains 

0^0 + a (8) 
For small distortion angle one gets 

cos2(0 + a) = (1 -2a 2 ) cos20-2a sin20 (9) 

The mean interaction potential per area is obtained from Eqs. [7| - ^ as 



u = 2^S 2 (l-2a 2 ) (10) 



Po 

In the case of a bend deformation, the distortion angle is on a molecular scale 

<90 

(id 

£ y represents the characteristic length of a bend deformation. One expects 
that £ y is the length of one elongated cell. The elastic distortion area density 
is then 

— 4 (l2 > 

which has to be compared with the phenomenological derived bend distortion 
energy 

^ 33 = 2 h3 ( iy ) (13) 



The splay elastic constant is then 



P 2 



k z , = U r -{SQ 2 (14) 
Po 

A similar expression can be derived for the splay elastic constant, k\\. One 
expects that the elastic coefficient increases with (i) increasing order pa- 
rameter, (ii) increasing cell density and (iii) increasing cell dimension. A 



10 



factor rs 500 is expected between the value at the nematic-isotropic transi- 
tion (p = HO cells/mm 2 and S ~ 0.5) and that of a nematic liquid crystal 
film of dense packed cells (p ~ 10 po and S ~ 1). 

The ratio of the elastic constants is determined by the ratio of character- 
istic length of the splay and bend deformation. 



hi = f£x\ 2 
k 33 UJ 



(15) 



One expects that these characteristic lengths, £ x and £ y , are approximately 
a. Only small differences between £ x and £ y are expected since a cluster of 
cells will change their orientation if one tries to change the orientation of a 
single cell. 

The coefficient b transforms the cell specific physical units into man-made 
units. Up to now no orientational elastic measurements were performed and 
consequently the coefficient, b, is unknown. It can be estimated in the follow- 
ing way: The orientational elastic coefficient close to the nematic-isotropic 
transition is expected to be the interaction potential, U, divided by the cellu- 
lar length, a, and multiplied by the thickness, d, of the nematic film (U -d/a) 
as shown above. The orientational elastic coefficient close to the transition is 
predicted by Eq. [14] to 2 6p £ 2 (with S ~ 0.5 at the transition). This leads 
to 

b^-—^U (16) 

One gets b ~ 5 x 10~ 20 J if one uses U ~ 5 x 10~ 19 J/cell as estimated above 
and d ~ 20pm, a « 100pm, p — 110 cells/mm 2 . 



Spontaneous Splay 

The spontaneous splay coefficient, k\ is zero if the axis perpendicular to n 
is a symmetry axis. If the building blocks carry a preferred direction in the 
direction of the long axis (e.g. wedge-shaped building blocks) there are as 
many building blocks up (f) as there are down (j) then the spontaneous 
splay coefficient is zero. A (parallel) polar nematic liquid crystal with its 
(parallel) polar distribution function (TtH) is formed in case of a (parallel) 
polar cell-cell interaction. The (parallel) polar distribution and the (parallel) 
polar shape asymmetry of the building blocks leads to a natural splay of the 
structure and the spontaneous splay coefficient, k\ is expected to be unequal 
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to zero. If the (parallel) polar anisotropy of the building-blocks (e.g. wedge- 
shaped cell, (parallel) polar organized cell structure) is connected with a 
(parallel) polar distribution of the muscle proteins then the center of mass 
of an area element moves in the direction of the director field. In case of a 
point defect divv ^ 0. 

To proceed further one has to introduce the motor characteristics of the 
migrating cells. The excess force, F(c, v), of a migrating cell depends on the 
concentration, c, of the migration stimulating molecules in the extracellular 



space and on the speed, v, of the cell [IS 



FM = F — (1-— ) (17) 

^ T -lie u max 

The first term describes the maximum force, F , which is 40 nN for a granu- 
locyte. The second term describes the fraction of membrane-bound receptors 
loaded with the migration stimulating molecule. The maximum speed, v max , 
is 24 /xm/min for granulocytes. One expects movements as demonstrated in 
Fig. |TJb. The hydrodynamic equations for nematic liquid crystals formed by 
migrating cells are not yet developed. 

Spontaneous Bend 

The spontaneous bend coefficient, ks, is zero if the axis parallel to n is a 
symmetry axis. If the building blocks carry a preferred direction perpendic- 
ular to the direction of the long axis (e.g. banana-shaped building blocks or 
elongated cells with the leading front on one side (e.g. keratinocytes)), there 
are as many building blocks right (— >) as there are left (<— ) then the spon- 
taneous bend coefficient is zero. A (transverse) polar nematic liquid crystal 
with its (transverse) polar distribution function (TtH) is formed in case of 
a (transverse) polar cell-cell interaction. The (transverse) polar distribution 
and the (transverse) polar shape asymmetry of the building blocks leads to 
a natural bend of the structure and the spontaneous bend coefficient, k^ is 
expected to be unequal to zero. Human keratinocytes (skin cells) plated on a 
surface migrate and cluster in large islands. The elongated keratinocytes in 
the cluster show the typical paving stone appearance |L9|| . If the (transverse) 



polar anisotropy of the building-blocks (e.g. banana-shaped cell, (transverse) 
polar organized cell structure) is connected with a (transverse) polar distribu- 
tion of the muscle proteins then the center of mass of an area element moves 
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perpendicular to the director field. In case of a point defect one expects 
curl z v 7^ 0. 



Disclinations 

The basic assumption of the used continuum theory concerns the smooth- 
ness of the director field, n(r). In practice, however, textures are observed 
originating from singularities in the orientational field. These discontinuities 
in the orientation are called disclinations. The energy of the disclinations 
themselves is unknown and the energy of the surrounding distorted director 
field can be calculated by making use of the elastic continuum theory. 
The director is given by 

n(r) — [cos<E»(a;, y),sin$(x, y)] (18) 

where <3>(x, y) is the angle between the director and the x-axis of a fixed coor- 
dinate system. In order to avoid complicated mathematics the one-constant 
approximation is used k = ku = k 33 . 

g = ^k(divn) 2 (19) 



The cylindrical coordinates (r, ip) are used to describe the singularity. 

(20) 



g(r,ip) = ^k 



'<9$\ 2 1 /<9$\ 2 



dr I r 2 \dtf) 



It is evident that this expression diverges at r = 0. The cause is a disconti- 
nuity in orientation at r = 0. 

The functional dependence of $ on r and ip is determined by requiring 
that the distortion energy must be minimum with respect to variations of $. 
This means that $ must satisfy the Euler-Lagrange equation. 

<9 2 $ 19$ 1 9 2 $_ o 
dr 2 r dr r 2 dip 2 

The trivial solution is $ = const, where the cells form an uniform oriented 
nematic liquid crystal. The general solution for the disclinations is 

$= m ^ + $ (22) 
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For an apolar nematic liquid crystal 2m is an integer and for a polar nematic 
liquid crystal m is an integer. Some examples are shown in Fig. |5|. 

The question whether the nematic liquid crystal formed by living and 
migrating cells has a polar or an apolar structure can be decided by investi- 
gating the types of defects, (i) Dense packed elongated melanocytes form an 
apolar nematic liquid crystal since the half- number disclination (m = —1/2) 
is found very often in the orientation pattern (Fig. |3|). However, this result is 



not specific for melanocytes. Other cell types like fibroblasts [18], osteoblasts 



(Fig. y), etc. show half- number disclination and, thus, form an apolar ne- 
matic liquid crystal, (ii) No half-numbered defects are found in the dense 
packed clusters of migrating granulocytes and keratinocytes. Granulocytes 
form a (parallel) polar nematic phase since this cell type like to create a 
defect with m — 1 and $o — (Fig. [I]). Keratinocytes form a (transverse) 
polar nematic phase since this cell type like to create a defect with m = 1 
and $ = 7r/2. 



Interacting Disclinations 

The deformation energy, W, of an isolated disclination can be obtained by 
integrating the elastic energy density in a circular area around the disclina- 
tion. 

W = W c + 7ckm 2 In (J^j (23) 

R, R c and W c are the radius of the circular integration area, the core radius 
of the disclination and the core energy of the disclination, respectively. The 
elastic energy, W — W c , depends on the type of disclination. The elastic 
energy of a ±0.5-disclination is a forth of that of a ±l-disclination. Therefore, 
the ±0.5-disclination should occur more frequently than the ±l-disclination. 
This prediction is in accordance with the observation. 

The deformation field of one disclination is altered if another disclination 
is introduced. The elastic interaction energy, Wu, is then [|l| 

(R \ 
IT) ^ 

where R12 is the distance between the two disclinations. As expected the 
interaction between disclination of opposite strength like mi = +0.5 and 
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m 2 = —0.5 is predicted to be attractive and repulsive for equal signs. The 
force between the two disclination is |l[ 

As expected the force is predicted to be large for small distances and small 
for large distances. 

The elastic interaction between two disclinations is orientation depen- 
dent: Two possible orientations in the case of a —0.5 disclination and a +0.5 
disclination where the interaction energy is minimum, are shown in Fig. [?]. 
The elastic interaction energy of the two orientations is the same in the case 
of equal elastic constants (k = kn = k 33 ) is used. But, in the general case 
(kn 7^ /C33) the elastic interaction energy is minimum for one of the orienta- 
tions |TJ|. The orientation as shown in Fig. |^a is energetically preferred in 
the case of kn > k 33 since a lot of bend deformation is between the disclina- 
tions. The orientation as shown in Fig. 0b is energetically preferred in the 
case of fcn < A; 33 since a lot of splay deformation is between the disclinations. 
The splay elastic constant, kn, is smaller than the bend elastic constant, k 33 , 
for an apolar nematic liquid crystal formed by melanocytes since in the case 
of a pair of ±0.5 disclination the orientation as shown in Fig. [?|b is prefer- 
entially formed. But up to now, we are not able to determine the relative 
value of the anisotropy, 2(kn — A; 33 ) / (fcn + ^33), of the elastic constants. 



Core of a Disclination 

The nature of the core of nematic liquid crystals formed by anisotropic 
molecules is not known. However, the core of a nematic liquid crystal formed 
by migrating cells can be observed directly in the light microscope. Different 
types of cores are observed as shown in Fig. ||: (i) The core of a disclination 
is an area free of cells. This type of core is frequently observed for very 
elongated cells like melanocytes (Fig. ||a). (ii) The core of a disclination is 
an area with isotropically distributed cell as shown in Fig. [5]b. This type of 
core is frequently observed for smaller cells like fibroblasts and osteoblasts, 
(iii) The core of a disclination is occupied by a star-shaped cell as shown in 
Fig. |8]c. The cells which form the nematic liquid crystal are in an elongated 
bipolar state. 

The defect in the orientation field attracts disturbances of the orientation 
field. One observes that other cell types or the same cell type with different 
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shapes are trapped in the core of a disclination. For example, melanocytes 
have usual two opposing dendrites. But some cells have three dendrites and 
these cells disturb the nematic phase and are, therefore, attracted by the 
orientational disclination (Fig. |8]c). 

In nematic liquid crystals formed by elongated molecules it is assumed 
that the molecules in the core of a disclination are in the isotropic state. This 
assumption can be verified for a nematic liquid crystal formed by elongated 
migrating cells. There exist a sharp boundary between the nematic and the 
isotropic phase. The nematic phase is destroyed if the elastic energy density 
is larger than the mean field. 

2 > 2b^- (26) 
Po 

The bend elastic constant divided by the calibration coefficient b can be 
expressed by Eq. [17]. One gets 

The maximal bending 30/ dy of a +1/2 disclination can be determined from 
photographic pictures. The characteristic length, £ y , of a bend deformation 
is expected to be in the order of the length of an elongated cell (£ y ~ a). The 
observed angular change by going from one cell to the next is approximately 
15 to 20 degrees. The angular change per cell as predicted by Eq. |2^ is 29 
degrees. The assumption that a strong bend deformation can destroy the 
nematic phase is verified since the predicted angular change is in accordance 
with the experimental observation. 



Outlook 

The practical importance of this study lies in its objective description of the 
liquid crystal phases formed by migrating and interacting cells. Two types 
of nematic liquid crystals, a polar and an apolar phase, are found to be 
formed by certain cell types in culture and such phases may be expected to 
be present in many other cell types. In addition, the objective description 
may be useful to define and elucidate cases of cellular dysfunction, and of 
altered cellular function induced by specific molecules like pharmacological 
active molecules, etc.. 
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The major objective of current and future research is to go beyond this 
phenomenological description of the thermodynamic phases formed by mi- 
grating and interacting cells. How the cellular machinery works on a molec- 
ular scale can be investigated by considering the technical data as a frame 
filled with the biochemical events. The cell, regarded as a fluid self-organized 
molecular machine, represents a new and important field in the triangle be- 
tween physics, chemistry and life sciences. The knowledge that one essential 
process in the cellular signal transduction chain is the supply of fresh recep- 
tors to the plasma membrane could lead to the application of new techniques. 
For instance the liquid crystal phase can be changed if special designed 
molecules are added which interact with the intracellular signal transduc- 
tion chain and alter cell-cell interaction either by the induced change in the 
cellular shape or by the transmitted chemical signal. 

The process by which migrating cells exchange information constitutes 
one of the most intriguing areas of life sciences and physics. An individual 
cell transmits signals and guides its migration and orientation according to 
the signals it receives, these messages being chemical, sterical or electrical 
in nature. The physical process of reducing a gas containing freely moving 
molecules to a liquid form is understood to a large extent and the mean 
behavior of a given molecule in the interaction field of the other moving 
molecules may be determined by Boltzmann statistics. Identical principle are 
used to reduce freely migrating cells to a condensed state. The mean behavior 
of a given cell in the interacting field of other cells can be described in terms 
of its automatic controller. The self-organization, involved in morphogenesis, 
organogenesis and wound healing, takes in a new light. 

The nematic liquid crystal formed by migrating cells opens new perspec- 
tives in the field of cellular tissue engineering. The two dimensional nematic 
liquid crystal can be used as a template. For example: (i) Cells of another 
type can migrate and orient on top of the liquid crystal to form more com- 
plex structures, (ii) Cells of another type can be attracted and localized in 
the orientational defects of the nematic liquid crystal. New types of liquid 
crystals are expected in analogy to blue phases where orientation defects are 
ordered in space, (iii) Complex structures formed by different types of cells 
can be altered if the phase of the template (cellular layer of the nematic 
liquid crystal) is altered by a chemical signal from a nematic liquid crystal 
phase to an isotropic phase. 
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Appendix 
Cell Preparation 

Melanocytes are one type of cells in the skin. They produce and distribute 
the color pigment melanin via dendrites. Cultured human melanocytes form 
the desired bipolar mode with two opposing dendrites when they are exposed 



to melanocytes medium. Details are given in ref. [20|. Briefly, skin biopsies 
were obtained from healthy donors. The cells were cultured in TIC medium 
(Ham F10, 16 % serum, 85nM PMA, 0.1 mM IBMX, 2.5 nM cholera toxin). 
Cells designated for the experiment were seeded in 25 — cm 2 dishes (Nunclon, 
Nunc, Germany) and cultured until confluence (approx. 25-30 days). 

The cells, being of spherical shape in the suspension, come down onto 
the substrate and start to develop their dendrites. They stay in this shape 
and crawl around on the surface. The cell body of different melanocytes 
attract each other but the dendrites repel each other. The orientation of the 
dendrites occurs during a periodic retraction and elongation process in such 
a way that each dendrite has a minimum of contact with the surrounding 
cells. 

We used a precision motorized x-y stage on a CCD-camera (Hamamatsu) 
equipped microscope (Zeiss Axiovert) and scanned an area of several mm 2 
using a lOx magnification lens. The pictures were digitized with a Scion frame 
grabber card on a Macintosh PPC and stored on the hard disk. The public 
domain software NIH-IMAGE (developed at the US National Institute of 
Health) was used for this purpose, as well as, for the further image processing. 

Picture Evaluation 

Our goal is to find an algorithm which determines automatically the mean 
orientation of the elongated cells in a small part of the picture. This task was 
successfully solved by going into the Fourier space. First, a squared section 
(64X64 pixels) of a picture is chosen, as shown in Fig. |9]a. A two-dimensional 
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Fast Fourier Transformation implemented in NIH-IMAGE is applied to this 
squared section of the picture. The result is a gray scale picture as shown in 
Fig. |b. 

The Fourier transform arises as a diffuse cloud with a handle-shaped 
contour. The maximum momentum of this elongated cloud is determined 
and found to be parallel to the direction of the mean cell alignment. In order 
to evaluate the exact orientation we apply a threshold procedure to remove 
all pixels with a gray value below a certainthreshold value. The remaining 
pixels are weighted with their gray value, g(x,y), and the angle, 9, of the 
orientation axis is found by maximizing the "momentum of mass". More 



details to the procedure are given in [21 



M xx -M yy + MM xx -M yy f + AMl y 

taut) = - (28) 

2M xy 1 ; 

with 

M xx = Y,x 2 g(x,y)- (Y, x 9(^y)f (29) 

and M xy and M yy analogous. 

This procedure is repeated for the next squared section of the picture. 
The scan was performed with overlapping squared sections. The shift length 
was half of the square size. 



References 

[1] P. G. de Gennes and J. Prost . The Physics of Liquid Crystals, Claren- 
don Press, Oxford University Press, Oxford (1995) 

[2] B. Alberts, D. Bray, J. Lewis, M. Raff, K. Roberts, J. D. Watson. 
Molecular Biology of the Cell, Garland Pub. Inc. New York, London 
(1994) 

[3] H. Gruler anf A. de Boisfleury-Chevance, J. Phys. I France 4, 1085 
(1994) 

[4] H. Gruler, U. Dewald, and M. Eberhardt to be published 

[5] P. R. Fisher, R. Merkl, and G. Gerisch, J. Cell. Biol. 108, 973 (1989) 



19 



[6] M. Schienbein, K. Franke, H. Gruler, Phys. Rev. E, 49, 5492 - 5471 
(1994) 

[7] R. Kemkemer, V. Teichgraeber, S. Schrank-Kaufmann, D. Kaufmann, 
and H. Gruler to be published 

[8] C. W. Oseen, Trans. Faraday Soc. 29, 883 (1933) 

[9] F. C. Frank, Disc. Faraday Soc. 25, 19 (1958) 

[10] E. Evans, D. Needham. J. Phys. Chem. 91 4219 (1987) 

[11] M. Schienbein, H. Gruler, Phys. Rev. £,52, 4183 - 4197 (1995) 

[12] H. Gruler. In: Biologically Inspired Physics. L. Peliti (Ed.), Plenum 
Press, pp 217 - 227 (1991) 

[13] H. Risken, The Fokker-Planck Equation, Springer Verlag Heidelberg 
(1984) 

[14] H. Gruler, Z. Naturforsch. 30a, 230 (1975) 

[15] Vereycken, H. Gruler, C. Bucherer, C. Lacombe, J-P Lelievre. J. Phys. 
Ill France 5, 1469 (1995) 

[16] I. E. Dzyalonshinsky, Sov. Phys. JETP 31, 773 (1970) 

[17] G. Friedel, Ann. Physique 18, 273 (1922) 

[18] R. D. Goldman, G. Berg, A. Bushnell, C-M Chang, L. Dickerman, N. 
Hopkins, M. L. Miller R. Pollack and E. Wang. In Locomotion of Cells 
Elsevier, North Holland Amsterdam (1973) pp. 83 - 107 

[19] L. Hegemann, A. Wevers, B. Bonnekoh, and G. Mahrle. J. Dermatol. 
Sci. 3, 103 (1992) 

[20] J. Griesser, D. Kaufmann, D. Maier, R. Mailhammer, P. Kuehl and 
W. Krone.J. Invest. Dermatol. 108, 275 (1997) 

[21] J. C. Russ, The Image Processing Handbook , CRC Press Boca Raton, 
Ann Arbor, London, Tokio (1995) 



20 




Figure 1: Human leukocytes exposed to blood plasma are migrating on a 
glass surface. The gap between the two glass slides is m 20fim. The fringe at 
the front of the cells is the amoeboid motor which pulls the rest of the cell, a) 
Single cell migration is observed at normal calcium concentration (2.5 mM). 
70 pictures were digitized in a time period of 390 s. The trajectories of the 
single cells are shown, b) The single migrating cells start to form a cluster 
of dense packed cells if the calcium of the blood plasma is bound by EDTA 
(< 5/iM calcium). The direction of migrating granulocytes is marked by an 
arrow. The immobile monocytes are marked by a circle. The monocycte 
close to the center of the picture acts as nucleation center of the migrating 
granulocytes. The polar order of the migrating cells in the cluster is very high 
(0.82). The monocyte at the upper right corner starts to act as an additional 
nucleation center. 



Figures 
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Figure 2: Human melanocytes on a plastic surface. The bar represents 
100 fim: a) Melanozytes spread at low cell density. Small sized clusters with 
oriented cells are formed, b) At high cell density the cell density becomes 
more uniform and the cells are oriented over large distances. 




Figure 3: a) Human melanocytes exposed to a smooth glass surface do not 
form an uniform apolar nematic liquid crystal, b) The calculated mean cell 
orientation of a) is shown. 
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Figure 4: A splay deformation is shown on the top and a bend deformation 
on the bottom. On the left side the change of the director described by ft is 
shown, if the space coordinate is changed by an increment dx or dy. The flux 
line of the director is shown on the right side (spay on the top and bend on 
the bottom). 




Figure 5: A schematic drawing of the flux lines of the director field is 
shown for four possible disclinations: a) and b) shows the diclinations for 
m = ±0.5 and $ — and c) and d) the disclinations for m = 1 and $ — 
and 7r/2. 
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Figure 6: A disclination (m = —0.5) is shown for osteoblasts. 




Figure 7: Two possible orientations of a +0.5 -disclination and -0.5- 
disclination are shown. The left side is preferred for ku > k 33 where the 
bend deformation is preferred and the ride side for ku < A; 33 where the splay 
deformation is preferred. 




Figure 8: A m = —0.5 disclination is shown for melanocytes. The bar is 
100 fim. Different possibilites are shown: a) The core of the disclination is an 
area free of cells, b ) The core of the disclination is an area with isotropically 
distributed cell, c) The core of the disclination is occupied by a star-shaped 
cell. The cells which form the nematic liquid crystal are in an elongated 
bipolar state. 
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Figure 9: A small rectangular frame of 64 x 64 pixels was chosen in the real 
space (photographic picture on the left side). The 2d-Fourier transformed of 
this picture is shown on the right side. The mean cell orientation is calculated 
by means of the Fourier Transformed and is shown by the bar. 
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